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Abstract 

We study the d - dimensional Bose gas at finite temperature using the renor- 
malization group method. The flow - equations and the free energy have been 
obtained for dimension d, and the cases d < 2 and d = 2 have been analysed 
in the limit of low and high temperatures. The critical temperature, the co- 
herence length and the specific heat of a two dimensional Bose gas have been 
obtained using a solution for the coupling constant which does not present a 
singular behavior. 
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I. INTRODUCTION 



Renewed interest for two dimensional dilute Bose gas has been recently determined by 
the discovery of high - temperature superconductors. 

The Bose - Einstein condensation is known not to occur at finite temperature for the 
interacting boson system in one and two dimensions (the Mermin - Wagner theorem); the 
absence of a condensate does not necessary imply the absence of a phase transition (in the 
universality class as XY model) and it is expected to be in the superfluid state. 

In the last years the study of two dimensional {2D) Bose condensation was stimulated 
by new experimental results concerning scattering in dilute gases of alkaline atoms, spin 
polarized hydrogen recombining and adsorbed helium monoalloys Jl|. 

The theory of weekly interacting Bose gas have been developed first at T = by Bo- 
golibov and Lee, Huang and Yang |3j at finite temperature. Using the T = results 
Hohenberg and Martin [|J, Popov and Fadeev |J, Singh || and Cheung and Griffin |7j 
developed the field - theoretical methods for finite temperature. Using this method it was 
possible to calculate the observable physical quantities, but these calculations are quite 
difficult because of the influence of correlations induced by the condensed phase. 

The many - body theory developed by Popov |||9] is manly based on the t - matrix 
approximation, also applied for the 2D Bose gas by Schick We mention that the t 

- matrix for the scattering in the Bose system involves a sum over all ladder diagrams to 
infinite order in interaction, taking the repeated scattering of two particles. 

The purpose of this paper is the study of the critical region of 2D interacting Bose gas, 
which was first done in the famous paper of Fisher and Hohenberg [12], using t - matrix 
method and the Renormalization Group (RG) method. The problem was also studied by 
Kolomeisky and Straley Jl3| for a D dimensional Bose system at T = 0. The influence 
of a gauge - field on a superfluid phase transition in a two dimensional model at finite 



temperature has been studied by Ubbens et al |14j and for the case of three dimensional 
dilute Bose gas by Biglsma and Stoof B. 
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As we mentioned, the results of RG study from Ref. 112] have been used for the problem 



of the quantum phase transitions in the spin systems |16|], or proposed |T7|] to study the 
large scale behavior in the crossover problem from high temperature superconductivity. All 
these papers used the approximations from Ref. |E? |, which will be reconsidered especially 
concerning the solutions of the RG scaling equations. We will also recalculate the number 
of condensate bosons for d = 2 and the coherence length using the new result obtained by 
solving more accurate the equations for the coupling constants. 

The paper is organized as follows. In Sec. II we present the model and write down the 
RG scaling equations. In Sec. Ill we solve the equations in the low temperature region and 
calculate the number of condensate bosons and the coherence length £(T) for the general 
case of dimension d. The interesting case d = 2 will be treated in Sec. IV where we also 
calculate the temperature dependence of the specific heat. The high temperature behavior 
will be analyzed in Sec. V. Sec. VI contain the relevance of our results compared with the 
results from literature. 



II. MODEL AND SCALING EQUATIONS 

We consider the dilute Bose system in d dimensions at finite temperature T described 
by the action: 



C _ c( 2 ) , c( 4 ) 
b *ff ~ O eff + b eff 



(1) 



where 



c( 2 ) _ " V" 
fc 



1 



2 7.2 



K k 
2m 



UJr, 



(2) 



and 



fcl &4 

and the following notation have been used: 



4 fci fc 4 Vi=l / 



(3) 
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?-^' T ?/#" (4) 

u) n = limksT being the bosonic frequencies. 

In Eq. (0) fi is the chemical potential of the bosonic system, described by the scalar field 
4>(k) and T is an energy parameter which controls the strength of the quantum fluctuations, 
the classical limit being T = 0. The Eq. (^) represent the interaction between fluctuations, 
u {u > 0) being the coupling constant. 

The basic propagator of the bosonic system has the form: 

'h 2 k 2 \Ur, 



G(k, lo t , 



(5) 



2m T 

and the renormalization transformations are carried out by integrating over a momentum 
shell and summing over all frequencies. 

Regarding this procedure we mention that the summation over the bosonic frequencies 
(transformed in integral) has to be done for all frequencies and not on a shell as in the 
method applied in ||l8lJl9l for bosonic excitations. The physical argument is that in this case 
the presence of the chemical potential is important because we expect different contributions 
on the energy scale, which is in fact identical with a change of physical behavior in different 
temperatures domain. On the other hand, the limit cases fi = ±oo correspond to the 
condensate, respectively normal states. 

After integrating out the degrees of freedom in the momentum shell, we rescale the 
variables as: 



and the field operator as: 



k — — , uj n = - — , T = — (6) 

b b z b z w 



i'(^, n ') = r^(p^) (7) 



where z is the dynamic critical coefficient and we will introduce for b the parameterization 
b = e l (b > , / > 0). The scaling equations obtained in Ref. ]r2| , |T4| for T ^ and in [T3| 
for T = will be written as: 



^ = -(2-^(0 



dT{l) 
dl 



zT(l) 



dfi(l) 
dl 



2M0-^[/x(/),T(/), M (/),r(/)] 



^ = [4 - (d + z)]«(o - ^ d {8F p _ p [/i(0, r(0, «(0, r(0] 
+2F p _ tt [^(0,r(0,u(0,r(0]}u 2 (0 



(d + s)F(Z) + K d F f [fx(l),T(l),u(l), r(/)] 



dF(Q 
d/ 

where F(Z) is the free energy and = 7r d / 2 / (2 d - 1 r(d/2)). The constants 
and from Eqs. (JX0| - 0) are given by the expressions: 

F M = F>(Z),T(Z),r(Z)] 

*QA 2 r(z) 



2/i(Z) - 



exp 



F p _ a = F p _ a [/i(Z),T(/),r(/)] 

1 



2r(Z) ( 



fi 2 A 2 
2m 



coth 



2A2 



r(z) /^a 



2k B T(l) V 2m 



-MO 



F p _ p = F p _ p [//(Z),T(Z)] 
1 



4 ^)-h 2 [^( 



' fi 2 A 2 
2m 



/i(Z) 



^ = ^1/1(0,7(0^(0] 

= k B T{l) In < 1 — exp 



r(z) /rA 



2A2 



- k B T(l) In 



1 — exp 



k B T(l) V 2m 

r(0 1 



/x(0 



Following the Ref. JT^JT^] we will solve these equations for d < 2 (e 
case d = 2 which is of special interest. 
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III. BELOW TWO DIMENSIONS 



In order to solve the scaling equations (H) - (12) given in the previous section we start 



with the equation (11) for the coupling constant at T = 0. Introducing the notation: 

mK d k d - 2 Y 2 {l) 



C 



2fi 2 



(17) 



this equation becomes: 



du(l) 
dl 



eu(l) - Cu 2 {l) 



(18) 



with the initial condition u(l = 0) = Uq. The fixed point of Eq. (fl8|) has the simple form: 



u 



and the general solution: 



«(0 



c 



V" 



(19) 



1 - C 1 - s) e -p (- 



(20) 



This solution satisfies the conditions: 



u{l = 0) = Uq 



(21) 



u{l — > oo) = u* 



(22) 



and in order to perform the calculation of physical quantities in the critical region we will 
use the linearized form of Eq. (|20|) 



u 



u(l) ~ u* H (u — u*)e 

Uq 



(23) 



In the limit of low temperatures the chemical potential /i(Z) given by Eq. (0) becomes: 

ri u(l')e- 2V 



HQ) = -K d \ d e 21 I dl' 
Jo 



exp 



ft 2 r(*')A 2 c _ 2 ;/ _ Tjl'MV) 



_2mk B T(l') e 



2mk B T(l') 



- 1 



(24) 



From Eq. (^) and (^) we obtain: 
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T(l) = Te- {2 - d)l 



(25) 



T(l) = Te 21 (26) 

and using for fi(l) the lowest approximation = fie 21 ) we calculate the second term from 
the exponential of Eq. Q2~3| ) as: 

2^T « 1 (27) 
This inequality is satisfied even in the low temperatures domain, because T _1 is the energy 
scale of the quantum fluctuations and in this domain these are important. Using these 
considerations, Eq. (|4j) we get: 

MO^-W/V ff? (28) 

\2mk B T ) 1 

and in order to calculate the condensed density n and the coherence length £ we will use 
Eq. ( p3[) for u(l). The renormalization procedure will be stoped for I = /*, given by the 
equation: 



2m 

where a ^ 1. If we introduce the notations: 

h 2 A 2 T 



= -J^- (29) 



2mk B T 



6 



(30) 



we obtain from Eqs. (EST) and 



-21* 



6=1 + - (31) 



1 + -) ~ -/(/*) (32) 



where 



„i* -26J' 
J(D = / dl' ; ^ (33) 

V ; Jo exp {ae- 21 ') - 1 V ; 



This integral can be expressed as: 

1(1*) c 

where A = 2e and 



A>4 

2a b A 



2 In A 



/oo , 



In the limit of low temperatures, Eq. fl32|) gives: 



e- 21 ' = 2e 



2mk B TT 

n 2 K 2 



1+e 



ini 

2e 



Using these result we can calculate the density of the condensed bosons as: 

d d k 1 



n = e 



-dl* 



{2n) d 



exp 



k B T(l 



- 1 



The integral from Eq. (|3 



„d-l 



A = K d dk 

Jo exp 



k B T(l 



- 1 



will be transformed using the notation 



2 L2 



2mk B T(l*) k B T(l*) 



as 



* = 2^ 



2mk B T(l* 



1 (d/2) 



|xm| e 1 - 1 



where |x m | = fi(l*)/k B T(l*). In order to perform the integral given by Eq. (PET) 
e small (d close to ci = 2) and in this case 



A ~ 



2mk B T{l* 



d/2 



hi 



1 — exp 



k B T{l*) 



Using T(Z*) = Te 21 * and Eq. (p9|) we obtain for the condensed density 



2mk B T 



In 



2eln 



2c 



s 



The temperature dependence of the coherence length £(T) will be calculated using 

r 2 (T) = ^Kni (43) 

where //(/*) is given by Eq. (P5|). This gives 



and using Eq. (|36|) 



Z(T) = -^-mjK-^T-^) (45) 



2eln ^V /2 (2m^)V2 



2c 



if we take a ~ 1. 

This behavior appear as universal and is not surprising because as it was mentioned in 



Ref. [13| the first contribution to the calculation of chemical potential in the fixed point 



gives a first term as universal. 



IV. TWO DIMENSIONS 



The dilute Bose gas at T = was studied in Ref. |13| and is known as marginal. In this 
case we consider the system at finite temperature in the low temperatures approximation. 
The free energy and the specific heat will be also calculated for this system taking z — 2. 
The high temperature limit will be also analised. 

A. Low temperature limit 

1. Scaling equations 
In this case the scaling equations have the form: 



(47) 



du(l) 



Mi) 



mK 2 2 

o~U (1) 



A 2 K 2 u 2 (l) 



dl ' exp ( 

The Eqs. ([46]) - © have the solutions: 

r(i) = r 



h 2 A 2 



2mk B T(l) 



T(l) = Te 



21 



«(0 



Anh 2 1 



m Z + Z 



where T will be take as T = 1 and L has been calculated as: 



lo 

The solution of Eq. ( f4"§| ) has the form: 

m(0 = 



Anti 2 
mu 



-21' 



2m Jo V + I 

which will be written as: 

4A 2 h 2 2l 2mk R T { 1 



( h 2 A 2 -9i' 

eX P {2^T e 



MO 



2m 



— In 



— (1 + - 

2Z V 



h A 2 { 2l 

\ -i 



2 A 2 



1 — exp 



h z A 



2mksT 



In 



2A2 



1 — exp 



h z A 



-21 



2mkBT 



2mkBT 
h 2 A 2 



F(l) 



where 



F(l) 



21 dx 



In 



2A2 



1 — exp 



h z A 



2mk%T 



lo (x + 2l ) 2 

The renormalization procedure will be stoped at Z* defined also by: 

h 2 A 2 



Mr 



— a- 



2m 



but in this case /i(Z) is given by Eq. (|55|) . Following the same procedure we calculate 
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-2; 



. 4 



a la a 



2mknT 1 



h 2 A 2 In?-* 



'A- 



(57) 



4 2mk B T 



and if we introduce the effective temperature T = h A 2 a/8mkB Eq. (|5TD will be written 
as: 

, r _ (7(a) T 1 



4 T lnf 



(58) 



Using Eq. fl44[ ) we get for coherence length 



lnf 



T_ 

To 



1/2 



(59) 



The number of condensate bosons has been calculated using a similar relation with fl3T|) and 
we get: 

1 



2mfc B T 
n(i J = -r^— in 



which gives 



4 ^ 2 l-exp(-5^e-»' 



n(T) = l^ lnln lT 



(60) 



(61) 



5. TTie /ree energy and specific heat 

The free energy given by general Eq. ([I]) has the form: 

dF{l) 



dl 



4F(l) + C f (v(l),T(l)) 



with 



r(/) fh 2 A 2 



knT \ 2m 



-m(0 



C/ = K 2 A 2 k B T In jl - exp 
Using the substitution T = Te 2x we write the solution of Eq. ( |62"D as 

F(Z) = jf' e-^C/ (Te 2x ) 

where we approximate Cj as: 



(62) 



(63) 



(64) 
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C f = K 2 A 2 k B T\n 
The expression for F(l*) becomes: 
F(l*) 

where 

.4 



1 — exp 



A 2 r r 
—k B T / die" 21 

2,71 Jo 



k B T 



ln(l 



-MO 



(65) 



e T 



(66) 



r h z A 



2a2 



kn \ 2m 



- V 



In the low temperature limit Eq. (4.21) will be approximated as: 



A 2 A 



T C(a) 
T 4 

lnf 



2 T 



hi 



hi 



A 2 

47T 



111 



(67) 



From this equation we calculate the specific heat C V (T) = —Td 2 T/dT 2 and the dominant 
contribution in temperature has the form: 



C V {T) = C 



In" 



(68) 



where C = Co (A) which shows that the result is A - dependent. 

B. High temperature limit 

The high temperature domain is called "classical" domain, because it is dominated by 
the classical fluctuations of the 0(x, r) and is usually defined by z = 0. 

1. Scaling equations 



From the general equations (|S|) - (|TTJ) we write these equations in the high temperature 
limit as: 
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d[T(l)T-H 
dl 



\T(l)T-\l)] 



^ = m)-^FMi),ni)Hi) 



(69) 
(70) 
(71) 



where 



v(l) = k B Tu(l) 



(72) 



and Fp, F v have been obtained from Eqs. (13|) - (|T5|) in the limit T — > oo as: 

A 2 



W).r(01 



£,[//(/), T(Z)] 



fi 2 A 2 
2m 



Ai(i) 



A' 2 



?i 2 A 2 



(73) 



(74) 



2m MO 

In order to solve the Eqs. ( |6~9|) - (f7T[) we have to define Z, the value of / at which the flow 
enters in the classical regime defined by 

T(l) 



r(0 



> i 



(75) 



and introduce 



fx(l) = no , u(l) = u , v(l) = v 

In order to simplify the calculation we perform a simple transformation I' 
make the flow to start at I' = 0. 

The new scaling equations describing the classical regime will be: 

d[T(l')T- l (l')] 



(76) 



I — I which 



dl' 



2[T(l')T-\l')} 



d»(i>) _ i AMO 



di' 



(77) 



(78) 



13 



[ ) TTh* A 2 



dV 



The equation ( j79|) has been obtained neglecting //(/') in Eq. 
can be solved and we obtain the exact solution 



(79) 



In this form the equation 



v(l') 



2v 



(80) 



Bv Q + (2 - Bv ) exp (-21') 

where B = 5m 2 /irk 4 A 2 . 

We define /'* a value of /' for which we stop the scaling by a similar condition as in the 
low temperature case, 



v(h) = i 



(81) 



which gives 



exp (2/',) 



2-Bv 
Bv - 2v 



(82) 



and from this equation we calculate 



l\ ~ I In — 
2 v Q 



(83) 



where we used uT(l'*) = 1. 

The Eq. ( ff8| ) for the chemical potential can also be solved using for v(V) the expression 
([SO]) and we get: 



A*(0 = e 2 '' 
where /2q is given by: 



Mo 



^f-^lnfe- 2 ''+ ^° 



2-Bv 



(84) 



m , m 



Using this equation /i(/') is approximated as: 



ln(e~ 2 ''*+ Bd ° 



2-Bv 



(85) 



//(0 ~ e 21 



2rnvo ( / _ , 

i-2 V * / 
7m 



(86) 
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2. Density of bosons 



In order to calculate the bosonic density, we will use the general equations: 

dF 



n = — 



(87) 



where the free energy will be written as: 

F(l') = — f dx e~ 2x T{x) In { 1 - exp 

2lX Jo [ 
Using the relations: 



T(x) I h X 



2\2 



k B T{x) \ 2m 



n(x) 



(88) 



n 



lo dfj,(l') \ dl' J dfi 
and for «(/') the simple approximation 

where u(Z) ~ «, the Eq. (|89"D becomes: 



(89) 



(90) 



n 



dl' e 2 ^) -±- 



(91) 



From Eqs. 



and ( |9~1~D we calculate the general equation for bosonic density n as: 



n = l_ fi'* ^ 

T 2rrJ 1 exp {r(h)/k B T(h) - n(h)} } - 1 
This integral will be approximated in the high temperature domain as: 

1 2m 



(92) 



n 1 r 1 '* 2m 
T 2n Jo h 2 A 2 



and from Eq. ( jS3[) we calculate: 



ft 2 A 2 „ . 
1 + ~2^ ( l} + 



2tt ^A 2 



(93) 



n m . 1 
~ In — 



(94) 



T 7rft 2 mo 

which is the density of bosons in the classical regime. This result gives us the possibility 
to perform a matching with the low temperatures regime. Indeed, using uq ~ 1/lnTo for 
Tq -C 1 we reobtain the expression for n c ~ Tlnlnl/T the well-known result from Ref. ||. 
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V. DISCUSSIONS 



We reconsidered the problem of dilute Bose gas at finite temperature using the RG 
method in d = 2. This problem has been considered first in the well-known paper by Fisher 
and Hohenberg fll2 |, but their calculation may give divergencies in the integrals contained 



in the physical quantities calculated in the quantum regime, because of an approximation 
for the coupling constant which is taken of the form v(l) ~ 1/7. 

We solved the problem by using a linearization of the exact solution near the fixed point 
for d < 2 and taking the exact solution in d = 2. We also calculate the coherence length £(T) 
and the specific heat from the free energy As a common point of view with the classical 
Fisher and Hohenberg paper we mention the non - universal behavior of the system. The 
new features which appear can be concluded as follows for d < 2 : 

• the bosonic density n(T) present the specific T d l 2 dependence as well as ln(l/elnl/e) 
dependence; 

l "I 1 /2 

• the coherence length £(T) calculated also has 1/ Ta^+f^eln 1/e dependence; 

For the two dimensional case d = 2 we used the exact solution for the coupling constant 
u(l) ~ 1/(Z + Z ) and we showed that : 

• the density n(T) is non - universal and n(T) ~ Tlnl/T, which gives for the critical 
temperature the Fisher and Hohenberg result, lnlnl/na 2 , where A 2 = 1/a 2 , a being 
the potential range; 

• the coherence length £(T) is non - universal and depends of T as £(T) ~ |(lnT)/T| 1 / 2 ; 

• the specific heat C(T) is also non - universal and C V (T) ~ T/| lnT| 3 ; 

The RG scaling equation have been used also in the high temperature domain for d — 2. 
The main result of this treatment is the calculation of the bosonic density which can be 
matched with the result from the low temperature regime by choosing the coupling constant 
as the solution of the scaling equations from the low temperature regime. 
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